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Nucleon selfenergies and spectral functions are calculated at the saturation density of symmetric nu- 
clear matter at finite temperatures. In particular, the behaviour of these quantities at temperatures 
above and close to the critical temperature for the superfluid phase transition in nuclear matter is 
discussed. It is shown how the singularity in the thermodynamic T-matrix at the critical temperature 
for superfluidity (Thouless criterion) reflects in the selfenergy and correspondingly in the spectral 
function. The real part of the on-shell selfenergy (optical potential) shows an anomalous behaviour 
for momenta near the Fermi momentum and temperatures close to the critical temperature related 
to the pairing singularity in the imaginary part. For comparison the selfenergy derived from the 
K-matrix of Brueckner theory is also calculated. It is found, that there is no pairing singularity in 
the imaginary part of the selfenergy in this case, which is due to the neglect of hole-hole scattering 
in the K-matrix. From the selfenergy the spectral function and the occupation numbers for finite 
temperatures are calculated. 
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I. INTRODUCTION 



Heavy ion reactions at intermediate energies probe the nuclear equation of state in a broad temperature and density 
range. One tries to extract from the observables signals for cluster formation, multifragmentation or the liquid- 
gas phase transition. The theoretical description of such phenomena demands to go beyond the usual quasiparticle 
description for the equilibrium properties (equation of state) as well as for the non-equilibrium properties (BUU- 
simulations). The need to go beyond the quasiparticle approximation in the description of nuclear matter is further 
advocated by the electron scattering experiments from heavy nuclei. These experiments give a clear evidence that 
the one-nucleon spectral function shows pronounced deviations from mean field estimates, which are due to NN- 
correlations jp . The effect of these correlations can most accurately be studied for nuclear matter j|] . 

A systematic quantum statistical approach for the description of dense nuclear matter can be given using Green 
function theory. For the description of phenomena like the formation of clusters (bound states) in the medium or 
the onset of a superfluid phase one has to allow for a finite lifetime (damping) of the one particle states. Such a 
treatment can be based on the nucleon spectral function, which is defined in terms of the real and imaginary parts 
of the selfenergy. Having the nucleon spectral function at the disposal, one is able to determine the momentum 
distribution, the response function as well as the nuclear equation of state including correlations beyond the mean 
field level. 

Most of the existing calculations of the nucleon spectral function have been done for zero temperature nuclear matter 
at the saturation density 72q. In refs. it is shown that the quasiparticle description is only valid near the Fermi 
energy. In general, the spectral function has a considerable width indicating that the one-particle states are strongly 
damped at the saturation density. Baldo et al. H calculated the on- and off-shell properties of the nucleon mass 
operator and the nucleon spectral function. They stressed the importance to retain the full frequency dependence 
of the selfenergy for the calculation of the spectral function. Kohler jij compared the nucleon spectral function 
calculated within Green function theory with a calculation done in Brueckner theory. He pointed out the near 
equivalence between the two theories at zero temperature. Vonderfecht et al. || used Green function theory to 
calculate the nucleon selfenergy within the ladder approximation at zero temperature. They emphasized the need 
to correctly treat the pairing correlations contained in this approximation if hole-hole propagation is included in the 
kernel of the vertex function (thermodynamic T-matrix). Benhar et al. || calculated spectral functions for nuclear 
matter using the hyper-netted chain (HNC) method. The spectral function and the nucleon momentum distribution 
in nuclear matter were calculated by Benhar et al. Q at densities below the saturation density within orthogonal 
correlated basis functions theory (OCBF). As a result, they find that with decreasing density the discontinuity at the 
Fermi surface decreases. This means that the momentum distribution does not approach the non-interacting response 
with decreasing density. This was interpreted as being due to the attraction in the N-N interaction leading to the 
formation of "bound pairs" of nucleons. Van Neck et al. Q use correlated basis functions theory (CBF) with the 
Urbana vn interaction to calculate the nuclear matter spectral function and momentum distribution below no. In 
agreement with the previous authors they found that the momentum distribution does not reach the non-interacting 
one with decreasing density. De Jong et al. calculated the nucleon spectral function based on an extension of the 
relativistic Dirac-Brueckner scheme. 

First finite temperature calculations of the nucleon self energy were carried out by Grange et al. |n|, Ropke et al. 
p| (in connection with the e.o.s.) and Kohler |12|. The latter author found a pronounced temperature dependence 
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of the spectral function at the saturation density. Using the extended quasiparticle approximation for the nucleon 
spectral function Schmidt et al. ]T^| included the formation of two-nucleon correlations (in particular bound states) 
in the equation of state of nuclear matter. The finite temperature nucleon spectral function in the low-density region 
of nuclear matter was calculated by Aim et al. jl4j demonstrating the necessity to include the contribution of two- 
nucleon bound states to the one-nucleon spectral function at low density. These authors find that genuine bound state 
formation in nuclear matter is only possible at extremely low densities (n < 0.05 no, T = 10 MeV). With increasing 
density the bound states are dissolved in the medium (Mott effect). However, the formation of two-nucleon pairs 
in the continuum as well as their possible Bose condensation |L5| can take place at these higher densities as well. 
The necessity to include the deuteron singularity in the two-particle T-matrix into the description of nucleon-nucleus 
scattering a low energies within the folding model was demonstrated by Love et al. |l6f| . 

Within this paper we would like to discuss the nucleon selfenergy, spectral function and momentum distribution at the 
saturation density and at finite temperatures. In particular, we will concentrate on low temperatures near the critical 
temperature for the superfluid phase transition in symmetric nuclear matter, which has a maximum of ~ 5 — 6 MeV at 
densities below no . This phase transition has recently been discussed by some authors in particular with respect 
to the possible formation of a condensate of neutron-proton-pairs in the 3 Si — 3 Di-channel (l?]] and the relation to 
the Bose condensation of deuterons in low-density nuclear matter [^8|. Moreover, Baldo et al. jl9| suggest that this 
could be a possible mechanism of deuteron formation in heavy-ion reactions at intermediate energy. The onset of the 
superfluid phase is contained in the thermodynamic T-matrix in the ladder approximation fl5f| . It manifests itself as a 
pole defining the critical temperature for the thermodynamic T-matrix [ ^o| . It could be shown that already above the 
critical temperature the T-matrix shows a resonance-like behaviour which could be understood as a precursor effect of 
the superfluid phase transition Consequently, near the critical temperature the selfenergy and correspondingly 
the spectral function should reflect this resonance-like behaviour of the T-matrix. In this paper we will demonstrate 
how the T-matrix, the nucleon selfenergy and spectral function are changed when approaching the critical temperature 
for the superfluid phase transition from above. 

In chapter II we derive a set of selfconsistent expressions for the retarded nucleon selfenergy and the nucleon spectral 
function at finite temperature within the framework of Matsubara Green functions. In the following chapter III the 
differences to the usual Brueckner theory generalized to finite temperatures are pointed out. The chapter IV contains 
results for the T-matrix and the K-matrix and the corresponding selfenergies including the optical potential at finite 
T. Within chapter V the corresponding spectral functions and occupation numbers are presented. 

II. GREEN FUNCTION FORMALISM AND APPROXIMATIONS 

For the derivation of the selfenergy and the spectral function we use the Matsubara Green function technique as 
outlined in ref. |2^ |. We put h,ks = 1 throughout the paper. From the definition of the one-particle spectral function 
A(l,oj) = i[G(l,cj + iO) — G(1,uj — i0)) and from the retarded selfenergy E the spectral function according to the 
Dyson equation reads 

Art \ 2ImE(l,w) 

A(l i u)) = - 2 , (1) 

[c-!^-Re£(l,c)P + [Im£(l,c)P 

where 1 = {pi, <Ji,t{\ denotes momentum, spin and isospin of a single particle. 
The one-particle spectral function (|]) fulfills the sum rule 
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(2) 



It determines the macroscopic properties of the system such as the occupation n\ of single particle states given by 

duj 



ni(n,T) 



2tt 



(3) 



where f(u>) = {exp((w — (i)/T) + 1} 1 is the Fermi distribution. The equation of state of nuclear matter, where the 
nucleon density n is a function of the chemical potential /i and the temperature T, reads 



(4) 



with f2 being the normalization volume. 

For the evaluation of the selfenergy X in eq. (Q) approximations have to be made. We start from a cluster decomposition 
of the selfenergy p3], which is given by the following set of diagrams 



m,zv) 




irred 



(5) 



where z u is the Fermionic Matsubara frequency. The cluster decomposition is appropriate for the consideration of 
n-particle correlations via the n-particle T-matrix T n . Within this paper we restrict to n — 2, i.e. to two-particle 
correlations represented by the first diagram in eq. (|5|). Within the ladder approximation the two-particle T-matrix 
is given by 



1 



T(121'2', z) = V(121'2') + - 2^ y(1234)G 2 (3456,z)T(561'2',z) , 

3,4,5,6 



(6) 



where the quantity Gi is defined as the product of two full one-particle Green functions given in spectral representation 
as 



G 2 (12l'2' )2 )= / p. [ ^lzMzM i(1)U)i(2)U0 , u , fe , 



(7) 



Using the spectral representation for the one particle Green function as well as for the T-matrix p3| the first diagram 
in (JsJ) yields 

dE [JV ) + 5 (£)]ImT 0X (1212, E) 



S(l, W + iO) = i^ J ^A(2,u>') (/(u/)F ex (1212) - J ^ 



E - w' - u - iO 



(8) 



where f(u>) is the Fermi distribution and g(E) — {exp((£' — 2/x)/T) — the Bose distribution function for the 

two-nucleon states. T cx and V ex denote matrix elements of the retarded T-matrix and of the potential including 
exchange. 



The ladder T-matrix equation (|6j,[7|) as well as the selfenergy (|8|) contain the one-particle spectral function (|l|) and 
thus form a selfconsistent set of equations. A solution of the complicated set of equations (0), (§,0) and (|J) can 
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be achieved by iteration. In the first step of iteration we replace the spectral function in eq. (^) by a quasiparticle 
spectral function 

A QP (l,u J ) = 27:S(u J ~e 1 ) , (9) 
where the quasiparticle energy ei is defined as 

2 

ei = |^+Re£(l,w) | w=ei . (10) 

Inserting the quasiparticle spectral function (|^) in eq. (Q) results in the quasiparticle approximation for the two- 
particle T-matrix (^). Within these approximations the following expressions for the imaginary and real part of the 
selfenergy (||) are obtained 

Im S(l, w) = ^ ^[/(e 2 ) + -9( £ 2 + ^ ImT ex (1212, e 2 + u) (11) 

2 



and 

d£ g(E)ImT ex (1212, E) 



ReE(l,o;) = i^^/( e2 )ReT ex (1212, e2 + a;)-P ^ ■ 



7T E — e 2 — lu 



(12) 



where P denotes the Cauchy principal value. T cx ist calculated in the quasiparticle approximation. Rather then (|l 
we use in the numerical calculation the real part from the dispersion relation (Kramers-Kronig-relation) in the form 

ReE(l,c) = £ HF (1) + P f ^L^E^l , (13 ) 

where E HF (1) denotes the Hartree-Fock shift. This was numerically checked to be equivalent with the explicit formula 
©• 

The bare nucleon-nucleon interaction was approximated by a separable ansatz 

rank 

V^ L '(p,p') = w L m (p)\ aio w L ^{p') . (14) 

i d=l 

The T-matrix then can be given algebraically as 

Tt L \p,p', p > E )=J2 w ^p^ - J °( p > m^K.kw^p') , (15) 

ijk 



t ( p f\ f d3p L \ t ( x <j^/(Pi2 + pW(P/2^-p) > 

MP,Eh = J L "~ (p) g _ e(p/2 +p) _ e{p/2 _ p) , (16) 



where < ... > denotes the usual angle averaging in the Pauli operator. Having the selfenergy at our disposal the 
spectral function follows from equation (|l]). The quasiparticle energies e\ as defined by © were determined in 
Hartree-Fock approximation 

£ HF (l) = i£/(e 2 )F ex (l,2,l,2). (17) 

2 
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With the T-matrix (|T|) one is able to calculate the critical temperature for superfluidity using the Thouless-criterion 
ppf . It has been demonstrated in ref. [^0) that this coincides with the critical temperature found in BCS-theory. For 
the separable ansatz ( p^ ) it reads 

det[l - J a {P = 0, E = 2/i, T = T c )]« = . (18) 



The Thouless-criterion for nuclear matter has already been evaluated in ref. 15 



III. SELFENERGY IN BRUECKNER THEORY 



Brueckner theory is only applicable to zero temperature systems. The near equivalence to Green function results has 
been demonstrated in this case as already pointed out above. A T > extension of Brueckner theory was formulated 
by Bloch and De Dominicis p4[ . An obvious difference between the two approaches (Brueckner and Green function) 
is that the former only includes particle ladders while the latter also includes hole ladders in defining the effective 
interaction. It is of some interest to numerically study the effect of this difference. We therefore define a Brueckner 
K-matrix as in eq. (|) with T replaced by K. The propagator (^) is however modified by the replacement 

1 - - /(w') — (1 - /( W ))(l - /(</)) . (19) 



The selfenergy (Vb) in Brueckner theory is given by [£5[ 

Vb(1, «) = i (/(ca)^™(1212, e 2 + w) 

2 

+ H V K^cx 1234, e 3 + £4 — : -<Wp2,p 3 + P 4 ■ (20) 

il w + e 2 - e 3 - e 4 - it] J 

For the discussion below, let us denote the two terms in eq. (^0|) by Vg and V^. The second term, Vg, is often 
referred to as the Brueckner second order rearrangement potential. At a temperature T = it is zero for states above 
the Fermi-surface and negative for states below, while the first term is positive for states above and zero below the 
Fermi-surface. As explained below, the retarded selfenergy in (finite-temperature) Brueckner theory, denoted by Sb, 
is given by 

ImS B (l,w)=ImV B 1 (l,w)-Imy B 2 (l,w) , (21) 



ReE B (l,w) = Rel/ B (l,c i ;) = Rel/ B 1 (l,w)-P / BV ' ' . (22) 

J IT LU — LU 

This quantity is to be compared to the Green function E defined in the previous section. 

While Eb is calculated numerically according to eqs. (p0|), (|2l|), and (p2|), rewriting these expressions to parallel those 
in the last section would facilitate comparisons between the two formalisms. To this end, one can, invoking unitarity 
on the imaginary part of the second term of eq. ( ^o|) , write Im Vb as p5| 

I m y B (l,w) = ^^[/(e2)-3(e2,w)]Im^ ex (1212, e2 + W ) , (23) 



G 



with 

SM = = f(e 2 )e-^ . (24) 

g(e 2 + uj) + 1 

Notice that at 2[i — e 2 + to where g has a singularity g is finite. It was argued in ref. |2?j that the Brueckner Vb 
should be identified with the chronological potential S c which is related to E (2^j2^] by 

ImE c (l,w) = tanh(i/3(w-^))ImE(l,w) . (25) 

In our notations here, this says, for T > 0, one can obtain IhiEb from 

B tanh(i/3(cj — /i)) ^ 
Now eqs. ( p^ ) and ( p4[ ) also give a simple relation between the two terms contributing to Iml^(l,cj): 

ImT/ B 2 (l,w) = -e-^-^lmVi(l,u;) . (27) 

It is then obvious that the division by the tanh function, in obtaining Im Eb from Im Vb, is equivalent to just switching 
the sign of lmV^(l,uj): 

ImS B (l, W ) = (l + e *-' l) )ImF B 1 (l lW ) , (28) 



as stated in eq. (|21[) above. 

An equivalent way of expressing the foregoing correspondence between Im Vb(1,u0 and Green function quantities is 
identifying 2z Im and 2iImV B 2 with, respectively, E < and E > , the nonequilibrium extensions of which govern the 
loss and gain collision terms in a transport equation |2§f| . The on-shell version of this correspondence was pointed out 
by Cugnon et al pij. 



IV. RESULTS FOR THE T-MATRIX (K-MATRIX) AND THE NUCLEON SELFENERGY 

In this exploratory calculation we used a rank one separable approximation (Yamaguchi potential p3|) as well as a 



rank two parametrization of Mongan |31|. The formfactors of the barely attractive Yamaguchi potential are of the 
following form 

w a {p) = 2 A ° 2 , (29) 



where the coupling constant and the effective range are given by 

An 



12.3178 (MeVfnT 1 )^ a = 1 S ! 



o 

14.6988 (MeVfm- 1 )^ a = 3 Si ' 
7 = 1.4488 fm" 1 . (30) 

The formfactors for the rank two Mongan potential, which contains in addition to a long-range attractive a short-range 
repulsive term, have the same form as given in eq. (p9|). The corresponding parametrization is given in ref. ]3l}]. For 
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numerical convenience we restricted only to S-waves (^So^Si) which give the dominant contribution to the T-matrix 
at low energies. 

As we are concerned within this paper with the modification of the selfenergy near the critical temperature for the 
superfluid phase transition in symmetric nuclear matter, we first study the onset of superfluidity in the temperature- 
density plane of nuclear matter. In Fig. [I] the critical temperature for superfluidity ( |l8| ) (a — 3 Si) is given as a 
function of the density. The solid curve refers to the Yamaguchi potential whereas the dashed curve is for the Mongan 
interaction. 




FIG. 1. The critical temperature for superfluidity in symmetric nuclear matter ( 3 Si-channel) according to eq. (|l^) as a 
function of density for the Yamaguchi and the Mongan potential. The saturation density no is indicated. 



Due to the repulsive component present in the Mongan interaction the critical temperature is reduced compared to 
the barely attractive Yamaguchi potential. These curves are in qualitative agreement with the calculations done in 
ref. |l5[ |. In the following calculations we fix the density to n = rig (dotted line in Fig. [j] and vary the temperature. 

The key quantity for the calculation of the selfenergy within the ladder-approximation is the thermodynamic T- 
matrix. In Fig. ^ the imaginary and the real part of the diagonal matrix elements of the thermodynamic T-matrix 
(triplett-channel) of the Green function theory ( |l5| ) are given as a function of the energy argument at fixed relative 
momentum p and total momentum P = as well as at a fixed chemical potential /i. The imaginary part shows a 
zero, which is located at lo = 2/i independent on the temperature. This is due to the fact that the imaginary part of 
the T-matrix (|l5|,|l6|) is proportional to 



l-2f(e(p))=g- 1 (2e(p))f(e(p))f(e(p)) , 



(31) 



which is obviously zero at 2e(p) = 2/t for any temperature. With decreasing temperature the slope of the imaginary 
part of the T-matrix at lo = 2/i increases; at temperatures near T c (T c — 4.02 MeV in our case) it has a characteristic 
principal value structure. The real part develops the corresponding peak at lo = 2/i while approaching the critical 
temperature from above. Finally at T — T c the real part of the T-matrix diverges in accordance with the Thouless 
criterion |goj . 

The results of Fig. |^ can be compared with the T — results by Ramos et al. |32| . They also find a pairing instability 
at T — and lo — 2ep which is most pronounced for small total momenta. The difference compared to our results lies 
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in the fact that although their imaginary part goes to zero at u> = 2ep as well, it does not change sign at this energy 
as it does in our calculations. 
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FIG. 2. Real and imaginary part of the diagonal elements of the thermodynamic T-matrix ( |l5| ) in the 3 Si-channel as a 
function of energy for temperatures far above and close to the critical temerature T c = 4.02 MeV and for zero total momentum 
at the saturation density. The relative momenta were set to p — 0.5 and 1.0 fm - . 



This is due to the fact that in ref. [B2| the Galitski-form of the two-particle propagator in the kernel of the T-matrix 
is used |22j whereas we use the Kadanoff-Baym-form |33|], which results in the retarded T-matrix instead of the 
chronological one used in ref. In the limit T — > both differ only in the sign of the imaginary part for ui > 2ep. 
Consequently both definitions differ also in the sign of the corresponding selfenergies at T = 0. At finite T the 
retarded selfenergy discussed so far and the chronological selfenergy as derived from the Feynman-Galitski T-matrix 
H are related by eq. @. 
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FIG. 3. Real and imaginary part of the diagonal elements of the thermodynamic K-matrix ( 3 Si) from Brueckner theory 
plotted for the same parameters as in Fig. |^. 



However, with this difference in mind we can conclude that the pairing instability at T = T c is in qualitative 
agreement with the corresponding instability at T — observed in ref. p2[. However, a proper treatment of the 
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pairing correlations below T c demands the inclusion of a finite gap in the single particle propagators. This has to be 
determined consistently from a combination of the BCS theory with the T-matrix approximation. Such a treatment 
known as quasiparticle-RPA |34j] has recently been applied to the one-dimensional Fermi gas |35| . 

To compare with the Brueckner theory the respective key quantity is the K-matrix. In Fig. ^ the K-matrix elements of 
the Brueckner theory are given for the same set of parameters. The imaginary part of the K-matrix does not change 
sign in contrast to the T-matrix elements (Fig. ||). This is due to the fact that in this case the Pauli-blocking (|l9| ) 
is positive in contrast to the Pauli-blocking term (|3l|), which changes sign at 2e(p) — 2/i. For low temperatures the 
imaginary part is effectively zero for energies below this particular energy. For energies above this value a pronounced 
maximum develops. This gives the corresponding maximum in the real part of the K-matrix. However, no critical 
temperature can be found where the K-matrix diverges as found for the T-matrix (Thouless-criterion). Thus, the 
different Pauli operator ( |l9| ) in the K-matrix leads to considerable deviations from the T-matrix case in particular in 
the limit of low temperatures. 

In order to calculate the spectral function it is necessary to evaluate the off-shell selfenergy, which itself has some 
interesting features. In Fig. |] the imaginary part of the retarded selfenergy (|ll|) is given as a function of u> at pi =0 
and n — n$. The results are plotted for various temperatures. 
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FIG. 4. Imaginary part of the nucleon selfenergy (Jllj) at the saturation density no as a function of energy for momentum 
pi = 0. The selfenergy is given for several temperatures above the critical temperature for superfluidity (T c = 4.02 MeV). 
The chemical potential fi = —23.5 MeV (T = 4.1 MeV) and the location of the pairing peak at luq = 2/x — ei = 35.6 MeV are 
indicated. 



At the highest temperature T = 20 MeV (long-dashed curve) the selfenergy is rather smooth. It gives a non-zero 
contribution for energies u> > —200 MeV, develops a maximum near to = —100 MeV and then starts decreasing slowly. 
In the T = 10 MeV case we find two additional extrema: a minimum at the the chemical potential oj = fi w —23 MeV 
and a second maximum at cj = 2/i — e(pi) « 35 MeV. Decreasing the temperature further this behaviour gets still 
more pronounced. At the minimum the value of the imaginary part is drastically reduced; however it still has a 
finite value and will reach zero only in the limit T — > (see discussion of Fig. ||) . The second maximum gets more 
pronounced as well. At a temperature T = 4.1 MeV, close to the critical temperature T c = 4.02 MeV, one observes a 
pronounced peak. 

This particular behaviour can be traced back to the behaviour of the T-matrix ( |i~5| ) at low temperatures. It has 
been demonstrated in ref. plfl that ImT has a zero at the particular energy value z = 2/i for pairs with zero total 
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momentum (compare also Fig. |]). This compensates for the Bose singularity in the imaginary part of the selfenergy 
( |i"T| ) rendering Im £ at this energy finite. However, in addition a critical temperature can be found such that the 
real part of the T-matrix becomes singular at the same energy (Fig. |J). As has been shown in Ref. Q] this critical 
temperature coincides with the one for the supcrfluid phase transition in agreement with the BCS theory. In fact, 
the singularity in the T-matrix is nothing but the wellknown Thouless criterion pp[ | for superfluidity If this critical 
temperature is reached the above mentioned compensation does no longer hold (see also ref. pl[|) and the singularity 
in the T-matrix leads to a corresponding singularity in the imaginary part of the selfenergy. The singularity is located 
at an energy loq = 2[i — e\. This is readily to be seen if one restricts to the pole part of the T-matrix with total 
momentum P = 0. Then the integration over p2 in ( pd| ) can be carried out directly yielding a (5-peak in Im S at the 
energy lvq. This shows that the singularity in the imaginary part of the selfenergy, which occurs for T — > T c is a direct 
consequence of the pole in the T-matrix at T = T c , indicating the onset of superfluidity. 

In Fig. ^|we continued the evaluation of the imaginary part of the selfenergy for temperatures below T c , disregarding 
for the moment the pairing instability. 
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FIG. 5. The same quantity as in Fig. ^]for temperatures below the critical temperature T c = 4.02 MeV. Again the chemical 
potential \i = -22.2 MeV (T = 1 MeV) is indicated. 



This allows us to demonstrate that for T — > indeed the value of the imaginary part of £ at the minimum (cu — 
|U = —23 MeV) approaches zero, in accordance with zero temperature calculations of various authors P]36|. This 
zero of the imaginary part of the selfenergy at u = \i is a wellknown property at zero temperature leading to the 
fact that particles at the Fermi surface have infinite lifetime, i.e. they are good quasiparticles. According to the 
Hugcnholtz-van-Hove theorem the chemical potential coincides with the binding energy per nucleon. The empirical 
value for this quantity at saturation is Eg /A = —16 MeV. The quadratic dependence of Im£ around u = [i at 
T = demonstrated in ref. |37| is also found in our numerical calculations. With increasing temperature a non- 
zero value for ImS is obtained at ui = im, however up to temperatures T < 20 MeV (see Fig. [|) a pronounced 
minimum is reminiscent of the zero-temperature property. At all temperatures T < T c the pairing instability shows 
up at wo = 2/i — ei = 35 MeV. This pairing instability indicates the breakdown of the T-matrix approximation 
in the vicinity of uq at temperatures below T c . A consistent treatment requires the inclusion of the BCS-gap for 
temperatures below T c . Such a complicated calculation has not yet been carried out. The T = 1 MeV curve in Fig. || 
is in qualitative agreement with the calculation of ref. except at energies around u> « 35 MeV. There we find the 
additional peak due to the pairing singularity in the T-matrix discussed above. 
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FIG. 6. The imaginary part of the nucleon selfenergy at the saturation density calculated in Brueckner theory (eqs. (|2l|)) 
for the same parameters as in Fig. H. 



In Fig. H the imaginary part of the selfenergy Im Sb ( pll ) is given as a function of energy for the same parameters as 
in Fig. ||. For T — 20 MeV Im Sb is in qualitative agreement with the corresponding curve in Fig. [|. With decreasing 
temperature the qualitative behaviour is similar to the Green function case given in Fig. |J except in the energy range 
around ojq . Whereas a strong singularity is seen in Fig. ^ for temperatures close to T c only a small maximum is found 
in the Brueckner case for the same temperature. As discussed above this singularity is due to the pairing singularity 
in the T-matrix and occurs at the critical temperature defined by the Thouless criterion. This behaviour does not 
show up in the Brueckner theory. 



60 




co [MeV] 

FIG. 7. Direct comparison of the imaginary part of the selfenergy in Brueckner theory (dashed curves) and Green function 
theory (solid curves) for temperatures T — 20 and T = 4.1 MeV. 

Fig. ^ gives a direct comparison between the imaginary parts of the selfenergy, calculated in the T-matrix and the 
Brueckner approximations, respectively. The upper curves show the case at T = 20 MeV. While for energies below 
uj = —100 MeV and above u) — 100 MeV the two curves almost coincide, the imaginary part of the Brueckner 
selfenergy lies below the one from the T-matrix calculation in the energy range between —100 and 100 MeV. An 
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analogous result has als been found in the zero-temperature case by Kohler ju. In the T — 4.1 MeV case (slightly 
above T c ) we observe a similar relation between the two approximations. However, at energies around ujq the Green 
function selfcnergy shows the pairing singularity, which is absent in the Brueckner calculation. 

In Fig. H the same quantity like in Fig. [i] is given using the Mongan potential instead of the Yamaguchi potential 
used throughout the rest of the paper. The purpose is to demonstrate the influence of the repulsive part of the 
nucleon-nucleon interaction on the selfenergy. 
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to [MeV] 

FIG. 8. The imaginary part of the nucleon selfenergy at the saturation density for the Mongan NN-interaction. The critical 
temperature in this case has a smaller value (T c = 1.58 MeV) than for the Yamaguchi-potential. Like in Fig. ^ the selfen- 
ergy is given for several temperatures above and close to the critical temperature for superfluidity. The chemical potential 
fi = 14.68 MeV and the energy loq — 82.46 MeV is indicated. The inset shows the region around ujq in detail. 



The repulsion leads to a lower critical temperature (T c — 1.58 MeV) as compared to the barely attractive case 
(T c = 4.02 MeV). The temperature dependence of ImE does not change qualitatively compared to Fig. ||. Thus, we 
suppose that the behaviour discussed above also holds for more realistic interactions, such as the Paris-potential |3^] 
as used e.g. in ref. pj]]. 
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FIG. 9. The real part of the off-shell nucleon selfenergy calculated from the imaginary part according to eq. ( |l3| ) as a function 
of energy for n = no and the same temperatures as in Fig. W. 
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In addition to the imaginary part the real part of the selfenergy contains important information, its on-shell part 
defining the optical potential for the nucleon in nuclear matter. In Fig. ^|the energy dependence of the real part of the 
selfenergy as calculated from eq. ([l3|) is shown for the same parameters as used in Fig. |[ Again there is a relatively 
smooth behaviour for the = 20 McV case. With decreasingtemperature a second minimum is found for energies below 
uiq. Please note that for T — 4.1 MeV (close to T c = 4.02 MeV) a principle- value-like behaviour around ui — ujq is 
found. This a direct consequence of the corresponding pairing peak in the imaginary part of the selfenergy at lu — ujq 
(see Fig. [|). 

In Fig. [l(] the real part of the on-shell self energy i.e. the real part of the optical potential, is given as a function 
of the momentum p. The upper curve shows a calculation using only the first term in eq. ( |l2| ) . This is a standard 
approximation often used to calculate the optical potential in nuclear matter [jl3 
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FIG. 10. The real part of the on-shell nucleon selfenergy as a function of the momentum p for n = no and different 
temperatures according to eq. (|l^). The upper figure was a calculation using only the first term of eq. ( |l^ ) denoted as ReS 1 . 
The lower plot shows the result of the full expression. The Fermi-momentum is indicated as pp. 



With decreasing temperature one observes a pronounced minimum near the Fermi energy, which has first been 
observed in ref. Jl3| and related to the inclusion of hole-hole scattering in the Pauli operator. However, using the 
full expression of eq. (^) for the evaluation of the optical potential the behaviour at low temperatures is changed. 
A particular structure is found around momenta p — pf which is enhanced with decreasing temperature. If one 
studies the behaviour of the real part in detail, one finds, that the 3 Si-channel of the T-matrix is responsible for this 
anomalous behaviour. We suppose, that the pairing peak in the imaginary part ImE(p,u) at temperatures T — ► T c , 
present also at finite momenta p, leads to a corresponding principal- value-like structure in the real part of the optical 
potential. Restricting to the pole part of the T-matrix this behaviour can be shown using the dispersion relation 
between the real and imaginary parts of the selfenergy. Consequently for the on-shell selfenergy ReS this leads to 
the wiggle at p = pf ■ 

Fig. [□] shows the same quantity calculated from the Brueckner theory. The upper plot denoted as Re Eg shows the 
temparature behaviour of the first order Brueckner term (eq. (|2^)). For the lowest temperature (T = 4.1 MeV) one 
observes a plateau-like behaviour around the Fermi momentum. Using the full expressions (lower plot) the repulsive 
second order contribution leads to a different behaviour at low temperatures, which is characterized by a strong 
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enhancement for low momenta and a corresponding minimum for momenta p > p-p. 
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FIG. 11. The real part of the on-shell nucleon selfenergy calculated with the Brueckner K-matrix. The parameters are the 
same as in Fig. [l(j. While the upper figure shows the results only with the first term of eq. ( p2| ) denoted as Re Eg the lower 
part shows the results of the full calculation. 



A particular behaviour of the optical potential for momenta around p = pp has also been found in refs. j^,[iT|]. In 
ref. JllJ a plateau-like behaviour was related to the behaviour of the effective mass at the Fermi surface. In |}9| it 
was shown that a non-monotonous behaviour around pp (anomaly) was entirely due to the strong attraction in the 
3 Si — 3 Di channel and that it was enhanced if the density was decreased below n — uq. According to our understanding 
the anomaly observed in ref. p9] is probably as well due to the pairing instability 143] discussed above. 
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FIG. 12. The real part of the on-shell nucleon selfenergy calculated in T-matrix and Brueckner approximations, respectively, 
for temperature T = 30 MeV. 



In Figs. |l2[]lj we give a direct comparison between the Green function theory and Brueckner theory with respect to 
the real part of the on-shell selfenergy for different temperatures. The upper curves correspond to the first term in 
eqs. (lit) and (p2p, respectively. As the form of these expressions coincides, the differences between Green function 
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and Brueckner theory in this case are entirely due to the different Pauli operators in the T-matrix and K-matrix, 
respectively (compare eq. ( |l9| ) . The lower curves show the result of the full expressions (|l2|) and (p^ ) . In Fig. [l^ we 
give the results for temperature T — 30 MeV. The differences between the two approaches are not very pronounced 
for this temperature, except for low momenta. With respect to ReS 1 the Green-function-curve is slightly enhanced 
compared to the Brueckner-curve for momenta p ~ p-p and slightly reduced below. This behaviour is reversed for the 
lower curves. 




p [fm '] 

FIG. 13. The same figure as Fig. |l| for temperature T = 10 MeV. 

In Fig. |l3| (T — 10 MeV) the differences between the two theories are much more pronounced. The different form 
of the Pauli operarator results in a non-monotonous behaviour for the Green function curve showing a pronounced 
minimum around p = pp. In contrast the Brueckner curve is monotonously decreasing. A behaviour like this has first 
been observed in r cf. Q. For the full expressions the differences between both theories are less pronounced showing 
an enhancement of the Green function curve with respect to the Brueckner curve for momenta p < pp. 




p [fm '] 

FIG. 14. The same figure as Fig. |l| for temperature T — 4.1 MeV. 
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In Fig. [TJI the temperature (T = 4.1 MeV) is close to the critical temperature T c = 4.02 MeV. One observes basically 
the same behaviour as in Fig. [l3[ although the differences are still more enhanced. 

In Fig. [l5| the imaginary part of the on-shell selfenergy (optical potential) is given as a function of momentum for 
various temperatures. The upper graph shows the contribution of the first term in eq. ([ll|) only. 
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FIG. 15. The imaginary part of the on-shell selfenergy as a function of p at no for the same temperatures as in Fig. |1C| , 
according to eq. ([H]). The upper figure was a calculation using only the first term of eq. ([H]) denoted as Im E 1 . The lower plot 
shows the result of the full expression. The Fermi-momentum is indicated as pf. 



One observes a pronounced temperature dependence for momenta below p = 2 fm _1 leading to a decrease with 
temperature towards a minimum near p = p-p. The lower graph displays the full contribution for the imaginary 
part (|Tl|). Again, one notes a strong temperature dependence below p — 2 fm _1 . For momenta around the Fermi 
momentum a pronounced minimum is exhibited with decreasing temperature. The value of Im £ at p = pf tends to 
zero for T — > 0. 

40 




FIG. 16. The imaginary part of the on-shell nucleon selfenergy calculated with the Brueckner K-matrix. The parameters are 
the same as in Fig. The upper figure shows the results only with the first term of eq. (|2^) denoted as ImSg. The lower 
part shows the results of the full calculation (|2l|). 



17 



In Fig. HI the imaginary part of the on-shell selfenergy calculated in Brueckner theory Im Eb is given as a function 
of momentum for the same parameters as in Fig. 15, The upper plot displays the first term of eq. (|2l|). Again, we 
observe a strong decrease with temperature for momenta below p = 2 fm _1 as in figure [l5|. However, due to the 
use of the different Pauli-operator in the K-matrix Im E B is basically zero for momenta p < pf in the limit of low 
temperatures. No minimum can be observed. In the lower graph the full contribution of Im Eb (|l]) is shown. One 
notes, that the temperature dependence for momenta around pf is in qualitative agreement with the Green function 
result (lower curve of Fig. |l5|). However, for momenta p < 1 fm _1 the temperature behaviour is reversed compared to 
the Green function case. Please note, that it is in the same momentum range, where one notes pronounced deviations 
in the real part of the optical potential at low temperatures (compare lower part of Fig. |l4j). The Brueckner results 
can be compared with the calculation of Im Vb(1, £i) in ref. [fTo)] . Taking into account the relation between ImEe and 
Im Vb (see sect. Ill) one notes the qualitative agreement of the two calculations. 

Summarizing, one observes pronounced differences between the optical potential at low temperatures calculated in 
Green-function-theory and Brueckner theory, respectively. These show up at momenta below the Fermi momen- 
tum. For the full expressions (lower curves) these differences are due to higher order terms in the Green-function- 
approximation, not included in the second order Brueckner calculation (see also ref. p5|). 



V. SPECTRAL FUNCTION AND MOMENTUM DISTRIBUTIONS 



From the off-shell selfenergy the nucleon spectral function is calculated using eq. ([[]). In Fig. [17] the nucleon spectral 
function is plotted as a function of energy ui for zero momentum p\ — at a density n = hq and for the same 
temperatures as given in Fig. ||. For T = 20 MeV one observes a quasiparticle peak at energies uj = —150 MeV and a 
background contribution extending up to energies to ~ 200 MeV. 
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FIG. 17. The nucleon spectral function at the saturation density as a function of energy for momentum pi = 0. The 
temperatures are the same as in Fig. W. 



With decreasing temperature the quasiparticle peak is slightly shifted towards higher energies and its width is reduced. 
This reduction is due to the fact that with decreasing temperature the imaginary part reaches zero at higher energies 
(compare Fig. |i|). In addition a second maximum forms at lower energies. This can be compared with the zero 
temperature results in ref. H for the spectral function at p = 0. Using Green function theory they also arrive at a 
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spectral function with two peaks of comparable size which are located at approximately the same energies as given in 
Fig. [r| in the T = 4.1 MeV case. 

The temperature dependence of the spectral function is not as drastic as one could expect from the change in the 
selfenergy with temperature (Fig. ||). Near the critical temperature the tail of the spectral function at higher energies 
shows additional smaller maxima, which result from the pronounced structures in the real part of the selfenergy. 
These in turn are due to the singular behaviour of the imaginary part at lo = 2/i — e\. 




to [MeV] 

FIG. 18. The nucleon spectral function at saturation density and T = 10 MeV in the energy-momentum plane. 

In Fig. |l8| the energy and momentum dependence of the spectral function is given at n — n Q and T — 10 MeV. One 
observes that the double-peak structure found at p — vanishes with increasing momentum. A single maximum 
remains for pi > 0.7 fm -1 which can be identified with the quasiparticle peak. The width of this peak is reduced at 
p = Pf due to the minimum in Im £ at the chemical potential. For larger momenta the peaks are broadened again 
until for very high momenta the width is reduced again. The latter behaviour is due to the fact that for very high 
momenta the influence of the medium represented by the selfenergy becomes negligible. Please note, that for our 
choise of the nucleon-nucleon-interaction there is no high-momentum tail of the spectral function, because it is barely 
attractive. 

In order to demonstrate the influence of correlations on the nucleon occupation numbers, the spectral function can 
be used to determine this quantity. In Fig. [it] the temperature dependence of the nucleon momentum distribution 
(occupation numbers) n(jp) ^ is given at a fixed density n = no- The correlated occupation number (full line) is 
compared to the corresponding Fermi distribution function (dashed line). In the T — 5 MeV case we observe a 
strong depletion for momenta below the Fermi momentum with a value of n(p = 0) = 0.78 compared to 1 for the 
non-interacting case. Above the Fermi surface we find a corresponding enhancement of the interacting occupation 
numbers compared to the non-interacting up to about p = 2 fru~ . For T = 10 MeV the depletion is less pronounced 
(n(0) = 0.82). This tendency towards the non-interacting occupation numbers is continued is the case of T = 30 MeV. 
With further increasing temperature the interacting response approaches the non-interacting one. Using the Brueckner 
K-matrix in ref. |45| ] the finite temperature occupation numbers are evaluated which are in reasonable aggreement 
with our results as well with the calculations of ref. (l(J . 

In Fig. ^0| the density dependence of n(p) at fixed T (10 MeV) is shown. For the sake of a better comparability 
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the curves are normalized to p/pf- When going to densities above no (n/no = 1.32) one observes a lower depletion 
{n(p = 0) = 0.85) compared to n(p = 0) = 0.82 at n = tiq. 
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FIG. 19. Occupation numbers at saturation density for different temperatures. The dashed curves represent the uncorrelated 
occupation numbers, i.e. the corresponding Fermi distribution. 



The same tendency has also been observed in ref. [|43j], although at a higher temperature (30 MeV). Going to lower 
densities (n/n = 0.61) we find the astonishing result that the depletion at lowp is further enhanced (n(p = 0) = 0.78). 




p/p F 

FIG. 20. Occupation numbers at T = 10 MeV for density values around the saturation density. Again the dashed curves 
represent the non-interacting case. 

On the other hand this corresponds to the zero temperature results of ref. indicating that the momentum 

distribution does not approach the non-interacting one in the limit of zero density. In ref. [^) the depletion at T = 
and p = stays rather constant for densities between no and no/2 at a value of n(p = 0) = 0.86. The authors of ref. 

interpret this result as being due to the attractive part of the nucleon-nucleon-interaction leading to the formation 
of bound pairs at low densities. 
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VI. SUMMARY AND CONCLUSIONS 



Using the Matsubara Green function approach, selfconsistent expressions for the nucleon selfenergy and the nucleon 
spectral function for nuclear matter at finite temperature were derived. The selfenergy and the nucleon spectral 
function at the saturation density were calculated in first iteration starting from the quasiparticle spectral function. 
The variation of these quantities with temperature was studied for temperatures close to the critical temperature T c 
for the superfluid phase transition in symmetric nuclear matter. We found that approaching the critical temperature 
from above a singularity develops in the imaginary part of the selfenergy. It was shown that this singularity is a 
direct consequence of a corresponding pole in the T-matrix at the energy E — 2/i, which indicates the onset of a 
superfluid phase at T = T c pc| . Thus, the modification of the selfenergy near T c can be understood as a precursor 
effect of the superfluid phase transition in nuclear matter. Another effect which was discussed as being related to the 
pairing instability in the 3 Si — 3 Di-channel is the occurrence of a wiggle around p = p-p in the real part of the on-shell 
selfenergy, especially pronounced when approaching T c from above. A similar effect although less pronounced was 
also found in the Brueckner calculation. 

The temperature dependence of the spectral function has been investigated for temperatures above the critical temper- 
ature. Despite the strong modification of the selfenergy there is no such drastic modification of the spectral function 
when approaching T c from above. This is consistent with the fact that below T c the condensate part of the T matrix 
is proportional to the square of the gap and consequently vanishes at the critical temperature. 

The momentum dependence of the spectral function shows considerable deviations from the quasiparticle behaviour 
at small momenta, whereas the quasiparticle picture holds approximately for momenta around p = pp as well as for 
large p. 

The occupation numbers were calculated from the spectral function at some finite temperature and density. It could 
be shown that with increasing temperature the non-interacting occupation number is approached. The depletion of 
the occupation numbers is enhanced with decreasing density at finite temperature. 

At the end we would like to mention some open questions related to our calculation of the nucleon spectral function: 
The first question is related to the model interaction we used in our exploratory calculation of the selfenergy and 
spectral function. For the energy and momentum range investigated in this paper the important features of the 
selfenergy obtained using the simple model interaction of Yamaguchi type were also found using a rank two Mongan 
interaction. It is supposed that these features remain of relevance also for more realistic potentials fl38fl . 

The second question is related to the problem of selfconsistency. In principle, the spectral functions have to be iterated 
until selfconsistency is reached. It remains to be seen, to what extent the features of the first iteration obtained in 
this calculation will also be found in a fully selfconsistent calculation. 

The third question concerns a consistent description of the system below T c , where the consistent inclusion of a finite 
gap is necessary for the evaluation of the spectral function. In principle our calculation is restricted to temperatures 
above the critical temperature for superfluidity. The Thouless criterion indicates the instability of the normal quasi- 
particle state with respect to the onset of superfluidity. A consistent treatment below T c has to be based e.g. on a 
BCS quasiparticle basis with a finite gap. Up to now such a calculation has not been carried out for nuclear matter. 
Instead, in most of the approaches at zero temperature the implications of the pairing singularity for the selfenergy 
were neglected. However, in a series of papers Dickhoff et al. || stressed the need to properly take into account the 
T matrix singularity, discussed above, which is present at temperatures below T c . 
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In conclusion we evaluated the nucleon self energy and spectral function for finite temperature. We compared the cal- 
culations within the Green function approach with a finite temperature generalization of the Brueckner theory. Special 
emphasis was put on the behaviour of these quantities near the critical temperature for the onset of superfluidity in 
nuclear matter. Within the Green function approach the pairing singularity in the T-matrix at the critical temper- 
ature generates a corresponding singularity in the imaginary part of the selfenergy. The non-monotonous behaviour 
(anomaly) of the real part of the optical potential for momenta p ~ pf could also be related to the pairing singularity. 
The spectral function at finite temperature shows a complex energy dependence, which cannot be generated from an 
energy-independent width. 

All the features discussed above cannot be incorperated into a simple quasiparticle description. Thus, the nucleon 
spectral function should be the appropriate quantity for the description of hot and dense nuclear matter. 
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